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S_i ■ Abstract 

On" 

We study the asymptotic behavior of a special smooth solution y{x,t) to the second 
member of the Painleve I hierarchy. This solution arises in random matrix theory and in 
the study of Hamiltonian perturbations of hyperbolic equations. The asymptotic behavior 
of y(x,t) if x — > ±00 (for fixed t) is known and relatively simple, but it turns out to be 
more subtle when x and t tend to infinity simultaneously. We distinguish a region of 
. algebraic asymptotic behavior and a region of elliptic asymptotic behavior, and we obtain 

rigorous asymptotics in both regions. We also discuss two critical transitional asymptotic 
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, <~| regimes. 

1 Introduction 

We study asymptotics for a solution to the second member in the Painleve I hierarchy, also 
^ ■ called the Pj equation, 

x = ty- Qy 3 + ^{y 2 x + 2yy xx ) + ^y xxxx J . (1.1) 

This is an ODE in the x-variable, but the equation depends also on the parameter t. Given 
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t, general solutions to (jl.ip have an infinite number of poles in the complex x-plane, but in 
[TJ, a special real solution to (jl.ip was constructed which has no poles for real values of x 
and t, and which has the asymptotics 



y{x,t) = ^z \x\ 1/3 + (D(\x\- 2 ), &sx->±oo, (1.2) 
for fixed t G R, where zq = zq(x, t) is the real solution of 

4 = -48sgn(x) + 24zo\x\~ 2/3 t. (1.3) 
It is remarkable but well-known that y(x,t) is also a solution to the KdV equation 

Vt + VVx + j^Vxxx = 0. (1.4) 

In the case t = 0, the existence and uniqueness of a real solution yo(x) with asymptotics 
(|1.2p were proved by Moore in |25j . Thus any real pole-free solution y(x,t) to (jl.ip with 
asymptotics given by (jl.2p solves the Cauchy problem for the KdV equation with initial data 
y(x, 0) = yo(x). Since the Cauchy problem for KdV is uniquely solvable with initial data 
yo (see [21]), y(x,t) can be characterized as the unique solution to the Cauchy problem for 
KdV with initial data yo- The Whitham equations corresponding to this KdV solution were 
already studied by Gurevich and Pitaevskii in [18] and by Potemin in |26j . 
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The equation (jl.lj) appeared for t = in a work by Brezin, Marinari, and Parisi [3J, where 
the authors gave physical arguments supporting the existence of a real pole-free solution to 
(jl.ip (for t = 0) with asymptotics given by 

~ ^F^xl 1 / 3 , as x -> ±oo, (1.5) 

see also [5J. The interest in this solution was renewed after Dubrovin [T2] conjectured the 
existence and uniqueness of a real solution to (jl.ip which is pole- free for real x, t £ M. The 
uniqueness part of this conjecture remains open until now. Another part of Dubrovins conjec- 
ture suggests that y(x,t) describes the universal asymptotics for Hamiltonian perturbations 
of hyperbolic equations near the point of gradient catastrophe for the unperturbed equation. 
This family of equations includes among others the KdV equation and KdV hierarchy, the 
de-focusing nonlinear Schrodinger equation, and the Camassa-Holm equation. In the partic- 
ular case of the KdV equation, the conjecture was proved in [2]. Recurrence coefficients for 
certain critical orthogonal polynomials (related to unitary random matrix ensembles) have a 
similar type of asymptotics involving y(x,t) [2j[8]. 

In the two above-mentioned applications the P 2 solution y(x,t) describes a singular tran- 
sition (when letting x and t vary) between a region of simple algebraic asymptotics and a 
region of more complicated oscillatory asymptotics involving the Jacobi elliptic ^-function. 
One may thus expect that y(x, t) itself also exhibits two different types of asymptotics. How- 
ever in the asymptotics (jl.2p . there is no trace of elliptic or oscillatory asymptotic behavior. 
The question arises whether elliptic asymptotics can be observed when letting x and t tend 
to infinity simultaneously. In this paper we will prove that indeed, depending on the precise 
scaling of x and t, y(x, t) admits either an algebraic or an elliptic asymptotic expansion. 

The type of asymptotics we obtain depends on the value of s = x|t|~ 3//2 as indicated in 
Figure [TJ If we let x — > ±oo, t — > -co or if we let x — > ±oo, t — > +oo in such a way that 
s = xt~ 3 / 2 remains bounded away from the interval [— 2\/3, 2 ^p ], an algebraic asymptotic 
expansion similar to (jl.2p holds. The leading term of this expansion will be determined by 
the equation x = ty — \y 3 } which is equation (jl.ip if we ignore x-derivatives. If x — > ±oo, 

t — > +oo in such a way that s G (-2 v / 3,^) (and remaining bounded away from the 
endpoints), the asymptotic expansion for y involves a system of modulation equations and 
elliptic 8- functions. This type of behavior was already suggested in [Ml [26], see also the recent 
paper [16] , For s near — 2\/3 and near 2 ^^ , two different transitions in the asymptotics for 
y take place. 

Remark 1.1 Asymptotics for y(x, t) where both x and t tend to infinity in such a way that 
s is bounded can be seen as small dispersion asymptotics for y as a solution to the KdV 
equation. Indeed for e > 0, by (jl.4p . u(x,t) := e 2 / 7 y(e~ 6 / 7 x, e~ 4//7 i) solves the KdV equation 
normalized as follows, 

e 2 

u t + uu x + —u xxx = 0. (1.6) 

Theorem 1.2 (Algebraic region) Let y be the real pole-free solution to the P 2 equation 

defined before. Suppose that either t < 0, s £ M or t > 0, s£R \ [— 2\/3, 2 ^p ]. Let zo(s) be 
the real zero (which is unique under the above conditions on s and t) of the equation 

zl = 24sgn(t)z - 48s. (1.7) 

In the limit where t — > — oo, we have 

y(s\tf/ 2 ,t) = z 4^\t\ 1 ' 2 + 0(\tn, (1.8) 
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Figure 1: y(x,t) has elliptic asymptotics if (x,t) goes to infinity between the two curves and 
algebraic asymptotics if (x, t) tends to infinity outside this region. 



uniformly for s£l. Ift—> +00 and s£l\ [— 2\/3 ; 27 5 L expansion $1.8]) holds as well, 
uniformly for s bounded away from the interval [— 2\/3, 
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Remark 1.3 If we only take into account the terms without x-derivatives in (jl.ip . we obtain 
x = ty — |y 3 . Changing va 
will be convenient later on. 



x = ty — |y 3 . Changing variables y 1— > 4f|t| 1//2 , this leads to (|1.7p . This change of variables 



For s € (— 2y/3, 27^ )1 the asymptotics for y(x,t) involve the third Jacobi elliptic ^-function, 
elliptic integrals, and a system of modulation equations given by 

(Pa + P2 + Pi? + 2{pl + P\ + Pl) = 120, (1.9) 
(& + P2 + Pi? ~ 4(/3f + /3 2 3 + Pf) = 360s, (1.10) 

1 

- &(£ - a) Vfc - tiVPi - £<% = 0, a = - + & + (1.11) 

where we are interested in solutions P$ < a < P2 < Pi- F° r the interested reader, we mention 
that the variables Pj = fa, a = t^a solve the elliptic Whitham equations appearing in the 
asymptotic theory for the KdV equation [TTJ, [271 E3 EH] in the particular case of initial data 

fa(x,0) = fa(x,0) = fa(x,0) = -(48^)1, 
see |26t [T7] . Those Whitham equations are given by 

ot dx 3 fa + a 3 

with 

f-\ r < (~r r ~2 P2 -h 

a{o-) = -p 1 + (p 1 -p 3 )——, a =2 r , 

and K(a), E(a) are the complete elliptic integrals of the first and second kind. Solvability 
of the system (|1.9p - (jl.lip was obtained by Potemin [26] for s 6 (-2^3, ^M), but can also 
be decuced from the solvability of the Whitham equations j27J [TTJ. The asymptotics for 
y(x, t) in the elliptic region are expressed in terms of the solution to this system. We obtain 
asymptotics of the same type as in |10| . although only rapidly decaying initial data for KdV 
were considered there. 
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Theorem 1.4 (Elliptic region) Let s G (-2\/3, 2 -^p) and let f3 3 < a < (3 2 < Pi solve the 
system of modulation equations il.9\) - [JTTl\) . Then y(x,t) has the following expansion if we 
let t — > +00, 

y(st^,t) = ^ 2 ^f^ + t 1 / 2 (/3 1 -/3 3 )|^ + ^(log^ (j^A +0{r^). (1.12) 

Here 9 is the third Jacobi 9-function 

+00 

0(z;t)= e 2wim2+ ™ 2 , (1.13) 



and 
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n = — VZ^te-^y/r^VPl^dZ, (1-15) 
15 ^/3 2 



K{cr) y /?i - #3 

and K(o~), E(a) are the complete elliptic integrals of the first and second kind. The expansion 
U.12\) holds uniformly for s £ (-2\/3 + 5, 2 4pp - 5) for any 5 > 0. 

If we take asymptotics for y(x, t) where t — > +00, x — > ±00 in such a way that s = xt~ 3 ^ 2 
tends to —2\/3 or 2 -^p sufficiently fast, neither the algebraic expansion nor the elliptic one 
is valid. In this case one can obtain two types of critical asymptotics. Similar transitional 
asymptotics have been proved for decaying negative smooth solutions to the KdV equation 
in [5j [6]. Because the proofs are very similar as for KdV and for conciseness of the paper, 
we will not prove those results here, but we will indicate below the procedure that can be 
followed in order to prove them. 

The first critical regime appears when s — > — 2v3, which corresponds to the left edge of 
the oscillatory zone in the (x, t)-plane, see Figure HJ The asymptotics involve the Hastings- 
McLeod solution [19J to the second Painleve equation. This is the unique solution to the 
Painleve II equation 

««(0 = fr(0 + 2g(fl 8 , (1.17) 

with asymptotics given by 

g(0~Ai(£), as £^+00, q (£) ~ 75 as £ _> _oo. (i.i 8 ) 
For £ G M, we have 

y(-2V3t 3 / 2 + c t 1/3 £, t) = 2V3t^ 2 lr&q(0 cos(i 7/4 u;) + 0(r 2 / 3 ), as t -> +00, 

(1.19) 

where Co = 5 1//3 , cj = 5 1/6 2 3 1/4 , and 

w = ^V5-3 3 / 4 -2-3 1 / 4 . 5 5 / 6 ^- 1 . 
21 
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In other words the leftmost oscillations can be modeled by a rapidly oscillating cosine and 
the amplitude develops proportional to the Hastings-McLeod solution. 



The second critical regime is the one where s — > 5 , i-e. the right edge of the oscillatory 
zone in the (x, i)-plane, see Figure HJ Here the asymptotic expansion for y involves a sum of 
sech 2 -terms: we have 



y 



9 7 00 

= -t-V^t 1 ' 2 + I>/l5t 1/2 X;sedi 2 (X fc ) + 0(t~ 5 / 4 In 2 t), (1.20) 
o o 

k=0 

r 3 

where C2 = 1 , and 
8-53 

X k = --(- + k - hit - In(VfcFfcfc) - (fc + ^) ln 7, (1.21) 

k 

h k = j 2 , 7 = 4v/2-5§ -7t .3"^. (1.22) 

7T* vfc! 

Recall that the KdV equation admits soliton solutions of the form asech 2 (6x — ct). The 
expansion (|1.20p can thus be seen as a superposition of solitons. The amplitude of the 
rightmost oscillations is of the same order as the leading order term of y outside the oscillatory 
region, and of the same order as the amplitude of the oscillations in the elliptic region. Because 
of the term with In t in (jl.2ip , every soliton is sharply localized near a positive half integer 
value of £. 

1.1 Riemann-Hilbert problem for P 2 

In order to prove our results, we will perform an asymptotic analysis of the Riemann-Hilbert 
(RH) problem associated to the real pole-free solution to the P 2 equation. For a more general 
description of Riemann-Hilbert problems for Painleve equations we refer to |15j . 

Consider the following RH problem for given complex parameters x and t, on a contour 
r = M~ U Uj=o e ^^^ consisting of eight straight rays orientated from left to right. 

RH problem for Y: 

(a) Y is analytic in C \ V. 

(b) Y satisfies the following jump relations on T, 

y + (C)=Y-(C)(j "A, for argC = ^ with j even, (1.23) 

F + (C)=F_(C)f 1 °\ for for arg ( = if- with j odd, (1.24) 

F + (C)=y-(C)(_° 1 J), forCGM". (1.25) 

(c) Y has an asymptotic expansion of the form 

y (0 = [l+ f>*0r, t)C*J C^Ne- e «^, as C -»• 00, (1.26) 
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Figure 2: The oriented contour T consisting of the four straight rays T\, T2, T3, and 



where 



iY = 7f(-i l) e_i ™ 3 ' ^(C;x,t) = I ^C 7/2 -^C 3/2 + < 1/2 , (1-27) 



and U3 



1 
-1 



This is the RH problem for a general solution to the Pj equation (jl.ip . The RH problem can 
only be solvable if the Stokes multipliers Sj satisfy the relation 

?) g t) a ?) Co t) (i ?) (0 t) a - (- . 9 ■ <-> 

For any set of Stokes multipliers sq, . . . , sg (independent of £, x, i) satisfying this condition, 
it was proved in [7] that the function 

y{x,t) = 2A 1M (x,t) - A 2 ljl2 (x,t), (1.29) 

with Ai given as in (jl.26p . is a solution to (jl.ip . However, y can have poles at certain isolated 
values of (x,t). Those points corresponds to the (isolated) values of x and t at which the 
RH problem for Y is not solvable. Since y solves (|l.ip . it also solves the KdV equation (jl.4p . 
Even in the more general situation where so, • • • , sq depend on £ (but not on x and t), y(x, t) 
still solves the KdV equation (locally near (x,t)) if the RH problem is solvable at (x,t). 

We are interested in one particular solution to the Pj equation. This solution corresponds 
to the case where 

si = s 2 = s 5 = s 6 = 0, s 3 = s 4 = s = 1, (1-30) 
which means that there are only jumps for Y on T = Uj =1 Tj, with 

ri = M + , r 2 = e^R + , r 3 = ifr, r 4 = e -^iR + , (1.31) 

see Figure [2j For this choice of Stokes multipliers, it was proved in [7] that the RH problem 
for Y is uniquely solvable for all real values of x and t, and that y defined by (jl.29p is the 
unique real pole- free solution to equation (11.11) with asymptotics given by (|1.2p . From now on, 
we refer to Y as the solution to the RH problem with Stokes multipliers given by (jl.30p . One 
can verify from the RH conditions that Y(Q has analytic continuations from each of the four 
sectors determined by Ti, . . . , T4 to the entire complex plane. We denote Yj, j = 1, . . . ,4 for 
the analytic continuation from the sector in between Tj and Tj mo< j 4+1. The jump conditions 
tell us that 

«=* -,(\ !). «=*(_° a J). !). y>-< 0- <-» 
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Outline 



In the next sections, we will obtain asymptotics for Y = Y(£; x, t) when x, t — > oo by applying 
the Deift/Zhou steepest descent method pXJ on the RH problem. Using (jl,29p . this will also 
lead to asymptotics for y(x,t). In Section [21 we will prove the algebraic expansion for y given 
in Theorem ll.2i In Section [3[ we will prove the elliptic asymptotics stated in Theorem 11.41 
In Section [U we will give a brief overview of the procedure that can be followed to prove the 
critical asymptotics (|1.19p and (ll.20p . In addition we will make some comments about special 
solutions to higher members of the Painleve I hierarchy and about unbounded solutions of 
the KdV equation. 



2 Algebraic region 

Here we generalize the asymptotic analysis 
[21[ I22j . We assume t ^ in what follows. 



done in [7] for t fixed, which was inspired by 



2.1 Construction of the (^-function 

In order to transform the RH problem for Y into a RH problem normalized at infinity and 
with 'sufficiently simple' jump conditions, we need to construct a so-called g-function. This 
function will be of the following form, 

g(() = ci(C - z ) 7 / 2 + c 2 (C - z Q f' 2 + c 3 (C - z Q f\ (2.1) 

where we take the branch of the roots which is analytic in C\ (— 00, Zq] and positive for £ > zq. 
Let us write 

s = x\t\- 3/2 . (2.2) 



If zq = zq(s) is the real root (which is unique for t < and for t > 0, s G E \ [— 2\/3, 2 27 5 ] ) 
of the third degree equation 

zl = 24sgn(t)z - 48s, (2.3) 

and if we let 

1 1 1 2 sgn(t) 

Cl= 105' C2 = 30 Z °' C3= 24 Z °-^' (2 - 4) 

it is straightforward to check that, with 6 given by (|1.27p . 

\t\ 7 / 4 g(C;s)=9(\t\ 1 / 2 C;x,t) + d 1 C 1/2 + 0(C 3/2 ), asC^oo, (2.5) 
where d\ does not depend on its explicit value can be calculated but is unimportant. 

Proposition 2.1 Suppose that either t < or t > and sGl \ [— 2\/3, Then 

g((;s)>0, for(>z , (2.6) 

Im<4(C;s)>0, for(<z . (2.7) 

Proof. The function (£ — zo)~ 3 / 2 g(( I ; s) is quadratic in £ and has no real zeros if c\ — 4ciC3 < 
0. Since c\ > 0, it is positive for large C and thus for all £ if c| — 4ciC3 < 0. Now by ([27 
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which is negative if either t is negative or t > and \zq\ > 4y|. If zq € (2\/2, 4y|), one 
verifies that both zeros are smaller than zq so that (j2.6|) still holds. In a similar way one 
shows that Im has no real zeros apart from zq if — 2IC1C3 < 0, which is true if 

f < and if t > 0, |«o| > 4\/3. Moreover in (— 4\/3, — 2\/2) both zeros lie to the right of zq 
and still holds. 

We can conclude that and (JUT]) hold for t < and for i > if z G R\ [-4-i/f, 4^3]. 
By (I23D this is equivalent to s £ R \ [-2^3, 2^]. □ 

By f)2. 1 j) and the above proposition, a straightforward complex analysis argument using the 
Cauchy-Riemann conditions leads to the following corollary. 



Corollary 2.2 Suppose that either t<0ort>0andseR \ [-2y/3, ^^]. For 9 > 
sufficiently small, we have 

Re 5 (C) < 0, for arg(C - Zq) = vr ± 6 , (2.9) 

Re 5 (C)>0, for arg(C - z ) = 0. (2.10) 

These inequalities enable us to transform the RH problem for Y to a RH problem for which 
the jumps decay exponentially fast to / as t — > oo, except on (—00,2:0) and in a small 
neighborhood of Zq. 

2.2 Normalization of the RH problem 

Fix < 9q < y- such that Corollary 12.21 holds. Let us define T as follows, 

r(c;M)= GiN i/4 1 



Vi(|t| 1 / 2 C;2;,t)el'l 7/4 ^ ;:c '*) CT3 , if < arg(C-2 ) < tt - 9 , 

Y 2 (\t\ 1 / 2 C;x,t)e\ t \ T/i ^,t)a^ if n _ 9q < arg (£ _ Zq ) < ^ 

y 3 (|t| 1/2 C;^,t)e l ' |7/49(C;a; ' t)CT3 , if vr < arg(C - 20) <^ + #o, 

Y 4 (|£| 1/2 C;M)e |i|7/4s(C;a; ' t)CT3 , if vr + ^o < arg(C - z ) < 2?r, 



(2.11) 



with c?i given by (|2.5p , and with the Yj 's the analytic extensions of Y as explained in Section 
1.11 Then T satisfies the following RH problem. 



2.2.1 RH problem for T 

(a) T is analytic in C \ S, where E = z + (RU e w±e °R + ). 

(b) T satisfies the following jump relations on S, 

T + (C)=T_(C)^ x J, for arg(C - Zq) = 0, (2.12) 

r + (C)=T_(C)f 2|t| r/4 (f . Xit) JV forarg(C-2 O )=vr±0 o , (2.13) 

^(0=^(0^ J), forarg(C-2 )=^. (2.14) 
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(c) T has the following asymptotic behavior as £ — > oo, 

T(o = (i + Bifo or 1 + o(r 2 )) |tr»T^Jv 5 (2-15) 

with 

Sl = ( d 4 + |t|- 1/2 Ai,n - di|t|- 1/4 A lil2 \t\- l ' 2 A hl2 - d x \t\-^\ ^ (216) 

where the values of the *-entries can be calculated but are unimportant for us, and 
with Ai as in (fL26|) . 

Using (jl.29p . we can recover y(x,t) from the identity 

y(x,t) = 2t 1 / 2 B hll -tBl l2 . (2.17) 

By Corollary 12.21 it follows that the jump matrices for T decay to the identity matrix when 
\t\ —> oo except on (—00, zq) and in a small fixed neighborhood of Zq. Ignoring a neighborhood 
of zq and ignoring exponential decay, we have jump conditions which can be solved explicitly. 

2.3 Outside parametrix 

We define the outside parametrix by 

P {co) (0 = \t\-^((-z )-^N, (2.18) 
which is analytic in C \ (—00, zq], and satisfies the jump condition 

4 oo) (C) = ^ oo) (C)(_° 1 I 

for £ G (— 00, zq) if N is given by fll .2TH . Note also that 

T(C) p(oo) (cr l = j + Cl + 0(c - 2); ag c ^ OQ) (2 1Q) 

with 

C 1 = B 1 - ^a 3 . (2.20) 
The leading order asymptotics for T and y(x, t) will be determined by the outside parametrix. 

2.4 Local parametrix near zq 

In order to obtain asymptotics for T uniformly for ( 6 C, we need to construct a local 
parametrix in a small fixed neighborhood U of zq. This local parametrix has been constructed 
in [71 Section 3.4] using the Airy function, and it solves the following RH problem. 

RH problem for P: 

(a) P is analytic in U\H. 

(b) P+(C) = P-(C) v t(C) f° r C ^ Sfl U, where vt is the jump matrix for T given by 
(p32l) - (f27[4l . 

(c) P(C)-P {oo) (C) _1 = f + 0(\t\- 3 / 2 ), as t -> -00 and if t -> +00 in such a way that s 
is bounded away from [— 2\/3, 2 "^F ], uniformly for ^ G 9?7. 

The local parametrix is needed for the rigor of the RH analysis, but its explicit expression 
in terms of the Airy function will not be needed; it does not contribute to the leading order 
asymptotics for y(x,t). 
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Figure 3: The jump contour T,r- 



2.5 Final transformation 

Define 

R(C;x,t) = < VS ' , S X 2.21 

\r(C;x,t)p(c)- 1 , forced 

Then one verifies that 

(a) R is analytic in C \ T, R , with T,r as given in Figure O 

(b) for C G Sfl, we have R + (() = R-(()v R ((), where v R (() = I + C(|t|~ 3/2 ) in the limit 
where t — > — oo or where t — > +oo and s remains bounded away from [— 2y/3, 
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(c) as C ^ oo, we have R(() = 1 + C(C -1 )- 

This is a small-norm RH problem which can be solved by a series expansion [9], and it follows 
from this standard procedure that 

R(C;x,t) = I + 0(\t\- 3 / 2 ), (2.22) 

uniformly in £ as t — > — oo and also as t — > +oo with s bounded away from the interval 
[— 2\/3, 2 "^p ]- Furthermore as ( -> oo we can expand R: 

R(() = I + ?± + 0((- 2 ), asC^oo, (2.23) 

which implies that 

Ri{x,t) = 0(|t|- 3 / 2 ). (2.24) 
Thus by ([2T20]) we have 



f^+o(itr 3/2 

which results in the asymptotics (jl.8p for y by (|2.17p . 



Bl = ^a 3 + 0(\t\- 3 / 2 ), (2.25) 
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3 Elliptic region 



In this section we consider the case where t > and s G (-2^3, ^P>). The construction 
done in the previous section fails in this case because Proposition 12.11 does not hold, which 
would result in exponentially growing jump matrices for T at certain parts of the contour. In 
order to prevent this, we are forced to modify the g- function and afterwards to open lenses 
not only along (— oo,zq) but also on an additional interval. 



3.1 Construction of the (^-function and modulation equations 

We search for a ^-function in the form 

9(0 = TjTr f (z - P 3 ) l/2 (z -a)(z- h) l /\z - P^ 2 dz, (3.1) 
analytic in C \ (— oo,/3i] and positive for ( > Pi, where 

a = ~(p 3 +p 2 + Pi), (3.2) 

with /?3 < a < @2 < Pi real and depending on s but not on As C — > oo, the ^-function can 
be expanded as follows 

g(C; x, t) = -U 7/2 + C1 C 3/2 + c 2 c 1/2 + d lt ~^c 1/2 + 0(C 3/2 ), (3-3) 

1U5 

with 

ci = - 3^((/3a + h + Pi? + ml + /3 2 + ?!)), (3.4) 

C2 = + fa + Pif - 4(/3| + p\ + ft)). (3.5) 

Now we want to choose P$, fa, and fa in such a way that 

t 7 ^g((;s) = 9(t 1 / 2 (;x,t) + diC 1/2 + 0(C 3/2 ), as C ^ oo. (3.6) 
This is true if 

(fa + fa + fa) 2 + 2{Pl + p 2 + P 2 ) = 120, (3.7) 

(fa + fa + faf ~ 4(/3f + Pi + Pi) = 360s. (3.8) 
Furthermore, in order to be able to do a steepest descent analysis, we also require 

rfo 

v/^ft^-^V^eVA^^^O, (3-9) 

ft 

or in other words g±(fa) = 9±(fa)- It will become clear later on why we need to impose 

the latter condition. We look for solutions to the modulation equations (13.7p - (13.9p for which 

fa < oi < fa < fa. For two special values of s, the system of equations can be solved 

easily. The first one corresponds to the confluent case where fa = fa = a. Then (|3.9p is 
automatically satisfied, and (I3.7p - (l3.8p are solved uniquely by 

fa = fa = ol = -V3, Pi = W3, s = -2V3. (3.10) 

In the case where fa = Pli we have the unique solution 

Ps = ~Vl5, « = -1^15, p 2 = p 1 = x/15, s = ^^- ( 3 - n ) 
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Figure 4: The jump contour E. 



Those confluent cases correspond exactly to the values of s which are at the border between 
the algebraic and the elliptic region. For s in between the two critical values —2y/3 and 2 ^p , 
the equations (|3,7p - (|3,9p are solvable [26]. Note that the different variable z = -4= was used 
in [26]. 

Before proceeding with the RH analysis, let us write down some useful properties of the 
g- function, relying on (13.ip and (13. 9|) : 



g+(Q + g-(Q = 0, 

g + (0 ~ 9-(0 = 2g + (/3 2 ) = 2g + ((3 3 ) 



-id, 



for Ce (-oo,/3 3 )U(/3 2 ,/3 1 ) 
forCG 



with 



9. 



15 J 02 



(3.12) 
(3.13) 



(3.14) 



3.2 Normalization of the RH problem and contour deformation 

We will now use the ^-function to normalize the RH problem at infinity in a suitable way. 
We define T in such a way that it has jumps on a lens-shaped contour as shown in Figured! 



T(C;x,t) :-- 



with d\ given by 



1 

dit 1 /* 1 



Y 2 (t 1 / 2 C,x,t)e tT/i 3(C^) 

y 3 (t 1/2 C;^,*)e' 7/45(C;x '' ) 



1T3 
<*3 

1 



Y 4 (t l / 2 (;x,t)e tV4 9^ 
. Then T satisfies the following RH problem. 



in region I, 
in region II, 
in region III, 
in region IV, 



(3.15) 



3.2.1 RH problem for T 

(a) T is analytic in C \ S, where E is as shown in Figure HI 
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(b) T satisfies the following jump relations on S, 

T + (C) = T_(C)(q x ). forC>/?i, 



:/ + ( O - '/'_ ( ,0 ( ^7/4^.^) 1 j , for C off the real axis, 



T+{Q = T.{q(^ 1 J), for C€ (-oo,/3 3 )U(/3 2 ,/3 1 ), 

T + (C)=T_(C)( ^ J, forCe(/3 3 ,/3 2 ). 

(c) T has asymptotics of the following form as £ — >• oo, 

T(C) = (/ + i)C -1 + 0(C" 2 )) H CT3 C ^ 3 iV. (3.16) 

Away from the branch points /?3 , /?2 ; /?i ? it is straightforward to verify that the off-diagonal 
entries in the jump matrices for T decay on £ \ ((— oo, /S3) U (/S2, /Si)) as t — > 00, if the lenses 
are chosen sufficiently close to the real axis. 

3.3 Outside parametrix 

Ignoring the small jumps and the branch points, we obtain the following RH problem. 

RH problem for p(°°) 

(a) is analytic in C \ (-00, 

(b) On (— 00, /Si), we have 



^(0 = ^(0 .! J' forCe(-oo,/3 3 )U(/3 2 ,/3i), (3.17) 

p(oo) (c) = p(oo) (c)e - ii7 /4c CT3; for C g (/ 3 3)A) . (3-18) 

(c) P (oo) behaves at 00 as 

p(°°)(C) = (J + DiC 1 + o(c 2 )) r^c^N. (3.19) 

Similar RH problems have been solved many times in the literature in terms of ^-functions 
and meromorphic differentials, see e.g. [10} [9j [H [23]. The only minor difference in our case 
is that 00 is a branch point. This RH problem can be solved explicitly using the Jacobi 
^-function and elliptic integrals. We will construct the parametrix explicitly and refer to pQ 
for the construction in more general settings. 

In addition to the RH conditions stated above, we need to construct the outside parametrix 
in such a way that 

p(°°)(C) = 0(|C-/?/ /4 ), as (^f3j,j = 1,2,3. (3.20) 

If this is not the case, the construction of local parametrices later on would fail. p(°°) has 
the form 

^-'(o = (^ l+ \)^<cr^( h f ,»,). (3,d 
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Here h, h are scalar functions we will determine below, and 

7(o=( (c -fvT - <3 - 22) 

If we construct h, h in such a way that they are analytic in C \ (—00, j3\] and that they satisfy 
the jump conditions 

MO = MC)e~ it7/4f \ for C G (ft, ft), (3.23) 

MO = MCK i7/4f \ for C G (ft, ft), (3-24) 

^±(0 = ^t(C), for CG (-00, ft) U (ft, ft), (3.25) 

it follows from (|3.2ip and (|3.22p that the jump conditions (|3.17|) - (|3.18p are satisfied. Moreover 
if h, h admit expansions of the form 

h(0 = 1 + hC 1/2 + hC 1 + 0(C 3/2 ), as C -> 00, (3.26) 

^(c) = i-/tir 1/2 +^r 1 +c(r 3/2 ) ) asc^oo, (3.27) 

and if we take d = hi, p(°°) has the expansion (|3.19|) with 

01,11 = h-h + P 1 + h 2 , Di, 12 = -hi- 1 ". (3.28) 

We will now construct h, h explicitly. Let us consider the third Jacobi ^-function 

+00 

0(z;t)= e 27rimz+7TiTm2 , (3.29) 



which is symmetric in z and has the periodicity properties 

6{z + 1;t)=Q{z;t), 9{z + r;r) = e- mr - 27riz 6(z;T). (3.30) 
Now we let h, h be of the form 

#(o)0MO + c) tr ,, fl(o) fl(c-tt(0) ,„„ n 



where 



"(0 = ^ / ^ (3.32) 
2C7 00 ^( z _ / 3 3 )( z - / 3 2 )(z-/3 1 )' 

C = l —;^^^^^^^=^^^^^^^= = = (3 33) 

7/33 V(z-ft)(z-ft)(z-ft) v^T 7 ^' 

• dz iK'(a) , d . 

r = i = — ^ / , 3.34 

7-oc CV(/3 3 -C)(/3 2 -C)(/3i-C) K(a) 

and is the complete elliptic integral of the second kind. We then have u+(C) = ~~ u -(() 

mod Z for C G (-oo, ft)U (ft, ft), and using (pT30|) we obtain (pT25l) . For ( £ (ft, ft) we have 
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u + (C) - u_(C) = T. Using the second property in (pOUl) . we obtain ([3^3]) - ([334") ) if c = 
with f2 as in (|3.14|) . Now we can calculate h\, h% in (|3.26p - (|3.27p . which leads to 

hi = h (? (0) _ i {c) ) ' (3 - 36) 

ft2 = ^U( c )-25(°) + ^ »-7< ¥ c »)- (3 ' 37) 

From the standard theory of ^-functions |13| it follows that 8(u(C)) has its only zero at 
C = /3 2 , so that P(°°) has no singularities other than the branch points. As £ — > /?3 , f3\ , h and 
/i is bounded and this leads to (1330]) . As C -> /%>, MC) and h{Q are of order 0(( - /3 2 ) _1/2 . 
Here we need to exploit the freedom to choose c\ in (|3,21[) in such a way that (|3.2U|) holds. The 
value of ci can be computed easily but is unimportant for us. This completes the construction 
of the parametrix. 

3.4 Local parametrices near (3 3 , (3 2 , and (5\ 

The local parametrices near the branch points can be constructed in the same way as in the 
algebraic case using the Airy function, we again refer to [9l [T5J E] for details. We do not need 
the precise form of the parametrices here, it is sufficient to have the existence of parametrices 
satisfying the conditions 

(a) P is analytic in a fixed neighborhood Uj of (3j, 

(b) P satisfies exactly the same jump conditions than T inside Uj, 

(c) For C € dUj, we have P(Q = P (oo) (C)(-f + 0(t-'^ 2 )) as t -> oo. 

(d) T(C)P -1 (C) is analytic at /3j. 

It should be noted that the construction of the local parametrix relies on fj3. 13f) and thus 
indirectly on ()3.9p . Condition (|3.20p is also crucial. 

3.5 Final transformation 

Define 

R((;x,t) = < VS ' , ' S X 3.38 

{T((;x,t)P(()-\ for (G ^ U C/ 2 U C/ 3 - 

Then after a similar argument as in the algebraic case, we have 

R{(;x,t) = I + 0(r 3/2 ), (3.39) 

uniformly in £ as t — > oo with s G (— 2\/3 + 5, 2 Y^ — 5), 5 > 0. As £ — > oo we have 

R(C) = I + Bl{X > t) - DliX > t) + 0(C% asC^oc, (3.40) 
which implies that 

Bi(x,t) = D 1 (x,t) + 0(r 3 / 2 ), (3.41) 
and by (pTTD and (^281) . 

= 2D lM t 1 l 2 -Di l2 t+0{t- 1 ' 2 ) = & ~ & + /3i t i/2 +2M i/2_^ t i/2 +0(t -l/ 2)- (3 42) 
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Using (|3.36|) - (|3.37|) . we obtain 



y ( x , t ) = P3 P^P^ t 1 ' 2 - t l ' 2 C- 2 {\og 9)"(0) + t l / 2 C-\\og 9)"(c) + 0(t- l l 2 ). (3.43) 
By standard manipulations for ^-functions and elliptic integrals, this leads to (|1.12p . 

4 Critical asymptotics for y and possible generalizations 

In this section, we indicate how the critical expansions (|1.19p and (jl.20p can be obtained. 
Afterwards we will make some remarks about asymptotics for the Brezin-Marinari-Parisi 
solutions to higher members of the Painleve I hierarchy and about asymptotics for certain 
unbounded KdV solutions. 

Painleve II asymptotics 

For s near we can proceed as in the algebraic case, see Section [21 with some mod- 

ifications. The first problem is that Proposition 12.11 does not hold for £ near 

— s/3. This 

leads to jumps on the lines zq + e 7r=te ' which are not uniformly close to the identity matrix. 
In order to overcome this, we need to close lenses again at the point —y/3, so that we have 
a contour as in Figure 0] but with (3% = fa. Then the jumps will converge to the identity 
matrix as t — > oo except in a small neighborhood of —y/3 < Zq. Near this point, we need to 
construct a local parametrix built out of ^-functions associated with the Hastings-McLeod 
solution to Painleve II [141 [20] . This construction is essentially the same as the one in [5] . 
The calculations that finally lead to (11.190 are rather tedious, and as they are similar as in 
[5], we do not think it is appropriate to include the details in this paper. 

Solitonic asymptotics 

For s near 2 ^F , we again proceed as in Section [21 but now Proposition 12.11 breaks down near 
vl5 > zq. There is no need to modify the jump contour for T here. It is however necessary 
to build a local parametrix near y/lE. This time the parametrix has to be constructed using 
Hermite polynomials, similarly as in [6]. The degree of the Hermite polynomials will depend 
on the value of £ in (|1.20j) . If £ is close to a half positive integer, a transition to Hermite 
polynomials of higher degree takes place, and this requires a modified local parametrix. A 
long calculation for which we refer to [6j leads to (|1.20p . 

Higher members of the Painleve I hierarchy 

The Painleve I hierarchy contains an infinite number of equations PJ™ of order 2m with 
m = 1,2... The equation for m = 1 is the Painleve I equation y xx = x + Qy 2 , and the 
equation for m = 2 is, up to a transformation x i— > ax,t i— > bt given by (II. 1|) . Brezin, 
Marinari, and Parisi [3] considered not only the case m = 2, but the general case m = 2k. 
They believed that for any k £ N, there is a real pole-free solution to the 2/c-th member 
of the hierarchy which has asymptotics of the form 

i 

Uk{x) ~ Tlcfc^l 2k + 1 , as x — > ±oo. (4-1) 

This conjecture was supported by Moore |25j when he considered the RH problem for the 
2/c-th member of the hierarchy. Although the general RH problem for the P 2k equation has 
4k + 2 Stokes multipliers, only three among them are non-zero for the special solution under 
consideration: the one corresponding to the positive real line and the ones corresponding to 
the two anti-Stokes lines closest to the negative real line, all three of them being equal to 
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1 with the orientation as in Figure [5J Adding 2k — 1 monodromy preserving time param- 
eters to the RH problem as in (|1.27p . where 9 now has a leading order term cC~^~ , one 
can study long time/space double scaling asymptotics as we did for k = 1. Following the 
general procedures of the Riemann-Hilbert analysis for Painleve equations [15] , one expects 
again regions of algebraic and elliptic asymptotic behavior, but in addition also regions of 
hyperelliptic behavior. The transitions between those regions would be interesting to study 
as well, and might lead to more general Painleve II hierarchy asymptotics, a more general 
form of solitonic asymptotics, and possibly also to critical asymptotics in terms of the special 
solutions to the equation with j < k. 

Unbounded solutions to the KdV equation 

We already mentioned that the RH problem for Pj generates solutions to the KdV equation 
by (jl.29j) . also in the case where the Stokes multipliers sq,...,sq depend on (. For generic 
choices of Sj((), y will have poles at certain values of x and t. However, for initial data 
satisfying y(x,0) = 0(|:r| 1//3 ) as x — > ±oo, the Cauchy problem for KdV is well-posed [21]. It 
would be interesting to see if such solutions can be generated by choosing appropriate Stokes 
multipliers Sj((), and if those solutions have asymptotic expansions similar to the ones for 
y(x,t). 

Acknowledgements 

The author is grateful to T. Grava and B. Dubrovin for useful discussions and comments. 
He is a Postdoctoral Fellow of the Fund for Scientific Research - Flanders (Belgium), and 
was also supported by Belgian Interuniversity Attraction Pole P06/02, by the ESF program 
MISGAM, and by ERC Advanced Grant FroMPDEs. 

References 

[1] P. Bleher, Lectures on random matrix models: the Riemann-Hilbert approach, 
larXiv:0801.1858l y2. 

[2] M.J. Bowick and E. Brezin, Universal scaling of the tail of the density of eigenvalues in 
random matrix models, Phys. Lett. B 268 (1991), no. 1, 21-28. 

[3] E. Brezin, E. Marinari, and G. Parisi, A non-perturbative ambiguity free solution of a 
string model, Phys. Lett. B 242, no. 1, (1990), 35-38. 

[4] T. Claeys and T. Grava, Universality of the break-up profile for the KdV equation in 
the small dispersion limit using the Riemann-Hilbert approach. Comm. Math. Phys. 286 
(2009), 979-1009. 

[5] T. Claeys and T. Grava, Painleve II asymptotics near the leading edge of the oscillatory 
zone for the Korteweg-de Vries equation in the small dispersion limit, Comm. Pure Appl. 
Math. 63 (2010), 203-232. 

[6] T. Claeys and T. Grava, Solitonic asymptotics for the Korteweg-de Vries equation in the 
small dispersion limit, arXiv:0911.5686. 

[7] T. Claeys and M. Vanlessen, The existence of a real pole-free solution of the fourth order 
analogue of the Painleve I equation, Nonlinearity 20 (2007), 1163-1184. 



17 



[8] T. Claeys and M. Vanlessen, Universality of a double scaling limit near singular edge 
points in random matrix models, Comm. Math. Phys. 273 (2007), 499-532. 

[9] P. Deift, T. Kriecherbauer, K.T-R McLaughlin, S. Venakides, and X. Zhou, Strong 
asymptotics of orthogonal polynomials with respect to exponential weights, Comm. Pure 
Appl. Math. 52 (1999), 1491-1552. 

[10] P. Deift, S. Venakides, and X. Zhou, New result in small dispersion KdV by an extension 
of the steepest descent method for Riemann-Hilbert problems. Internat. Math. Res. 
Notices 6 (1997), 285-299. 

[11] P. Deift and X. Zhou, A steepest descent method for oscillatory Riemann-Hilbert prob- 
lems. Asymptotics for the MKdV equation, Ann. Math. 137 (1993), no. 2, 295-368. 

[12] B. Dubrovin, On Hamiltonian perturbations of hyperbolic systems of conservation laws, 
II: Comm. Math. Phys. 267 (2006), no. 1, 117-139. 

[13] J. Fay, "Theta Functions on Riemann Surfaces", Lecture Notes in Mathematics 352, 
Springer Verlag (1970). 

[14] H. Flaschka and A.C. Newell, Monodromy and spectrum-preserving deformations I, 
Comm. Math. Phys. 76 (1980), no. 1, 65-116. 

[15] A.S. Fokas, A.R. Its, A. A. Kapaev, and V.Yu. Novokshenov, "Painleve transcendents: 
the Riemann-Hilbert approach", AMS Mathematical Surveys and Monographs 128 
(2006). 

[16] R. Garifullin, B. Suleimanov, and N. Tarkhanov, Phase Shift in the Whitham 
Zone for the Gurevich-Pitaevskii Special Solution of the Korteweg-de Vries Equation, 
larXiv:0912.48531 

[17] T. Grava and F.-R. Tian, The generation, propagation, and extinction of multiphases in 
the KdV zero-dispersion limit, Comm. Pure Appl. Math. 55 (2002), no. 12, 1569-1639. 

[18] A. V. Gurevich, L.P. Pitaevskii, Nonstationary structure of a collisionless shock wave 
Soviet Phys. JETP 38 (1974), no. 2, 291-297. 

[19] S.P. Hastings and J.B. McLeod, A boundary value problem associated with the second 
Painleve transcendent and the Korteweg-de Vries equation, Arch. Rational Mech. Anal. 
73 (1980), 31-51. 

[20] A.R. Its and V.Yu. Novokshenov, " The isomonodromic deformation method in the the- 
ory of Painleve equations", Lecture Notes in Math. 1191, Springer- Verlag, Berlin - New 
York (1986). 

[21] A. A. Kapaev, Weakly nonlinear solutions of the equation Pf, Journal of Mathematical 
Sciences 73, no. 4, (1995), 468-481. 

[22] A. A. Kapaev, Quasi-linear Stokes phenomenon for the Painleve first equation, J. Phys. 
A 37, no. 46, (2004), 11149-11167. 

[23] A.B.J. Kuijlaars and M.Y. Mo, The global parametrix in the Riemann-Hilbert steepest 
descent analysis for orthogonal polynomials, preprint arXiv:0909.5626 , 

[24] A. Menikoff, The existence of unbounded solutions of the Korteweg-de Vries equation, 
Comm. Pure Appl. Math. 25 (1972), 407-432. 



18 



[25] G. Moore, Geometry of the string equations, Comm. Math. Phys. 133, (1990), no. 2, 
261-304. 



[26] G.V. Potemin, Algebro-geometric construction of self-similar solutions of the Whitham 
equations, Russian Math. Surveys 43 (1988), 252-253. 

[27] F.-R. Tian, Oscillations of the zero dispersion limit of the Korteweg-de Vries equation, 
Comm. Pure Appl. Math. 46 (1993), 1093-1129. 

[28] S. Venakides, The Korteweg de Vries equations with small dispersion: higher order Lax- 
Levermore theory, Comm. Pure Appl. Math. 43 (1990), 335-361. 

[29] G.B. Whitham, Linear and nonlinear waves, J. Wiley, New York (1974). 
Tom Claeys 

Cite Scientifique - Laboratoire Painleve M2 
F-59655 Villeneuve d'Ascq, FRANCE 
E-mail : torn. claeys@math.univ-lillel . f r 



19 



